Noise calculations within the second-order von Neumann approach 
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We extend the second-order von Neumann approach within the generahzed master equation for- 
mahsm for quantum electronic transport to include the counting field. The resulting non-Markovian 
evolution equation for the reduced density matrix of the system resolved with respect to the num- 
ber of transported charges enables the evaluation of the noise and higher-order cumulants of the 
full counting statistics. We apply this formalism to an analytically solvable model of a single-level 
quantum dot coupled to highly biased leads with Lorentzian energy-dependent tunnel coupling and 
demonstrate that, although reproducing exactly the mean current, the resonant tunneling approx- 
imation is not exact for the noise and higher order cumulants. Even if it may fail in the regime of 
strongly non-Markovian dynamics, this approach generically improves results of lower-order and/or 
Markovian approaches. 
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Introduction. — Although electronic transport 
through interacting nanostructures has been a subject 
of intense study in the past decade, its theoretical de- 
scription in general situations remains a challenging task. 
Standard complementary theoretical approaches are the 
Landauer-Biittiker scattering formalism in the limit of 
negligible Coulomb correlations and the master equation 
for weak coupling to the leads, respectivelyi^i^ If both 
the single-particle coherence effects due to coupling to 
leads and many-body Coulomb correlations are impor- 
tant, highly complex phenomena can occur, such as, e.g., 
the Kondo effect at low temperatures.? Even well above 
the Kondo temperature the interplay between coherence 
and correlations may lead to nontrivial consequences?^ 
which are hard to describe theoretically. 

The second-order von Neumann approach (2vN) based 
on decoupling of equations of motion for the density 
matrix, ^'^ like its relative the resonant tunneling approxi- 
mation (RTA) based on a suitable infinite resummation of 
the perturbation series for the generalized master equa- 
tion (GME),^'^ is capable of capturing such effects, as 
has been confirmed also by the successful description of 
pertinent experiments.^'^ In Ref. 8 the line- widths of the 
nonlinear conductance for a quantum dot asymmetrically 
coupled to leads in the Coulomb blockade regime were 
measured for both voltage-bias polarities and fully ex- 
plained by the RTA, while in Refs. and[l3 a canyon of 
suppressed linear conductance through an InSb nanowire 
quantum dot was correctly reproduced by the 2vN as a 
function of the gate voltage and magnetic field. Further- 
more, in a series of worksj^iii the 2vN has been applied 
to the single resonant level as well as double quantum 
dot models, finding that corrections beyond the lowest- 
order sequential tunneling (cotunneling) were described 
very well and conjecturing that the method yields exact 
solutions for non-interacting models.^ This is a nontrivial 
feature as the method stems from the atomic limit (de- 
scribed by GMEs) complementary to the non-interacting 



case. 

Altogether, 2vN/RTA is a powerful and successful ap- 
proximation and it is, therefore, natural to ask whether 
it can be extended to the calculation of other transport 
quantities such as noise and full counting statistics. This 
is exactly the subject of this Brief Report — we extend 
the second-order von Neumann approach to the evalua- 
tion of the full counting statistics, and noise in particular, 
by incorporating the counting field in the theory and then 
test the performance of this generalization for the sim- 
plest possible model of a single-level quantum dot. We 
find that contrary to the stationary mean current, which 
is captured exactly by the 2vN, the noise and higher- 
order cumulants constituting the full counting statistics 
are not exact. Nevertheless, it does yield significant im- 
provement over simpler, more standard approximations. 

Generalized second-order von Neumann approach. — 
In our generalization of the second-order von Neumann 
scheme we use the equations of motion approach and 
closely follow the corresponding derivation without the 
counting field in Ref. [3. We only summarize the main 
steps and associated approximations and show the final 
results. A fully detailed derivation can be followed in 
Ref. [l^. We consider a generic transport Hamiltonian4 
consisting of three parts, H = i/systcm + ^^icads + Ht, 
with the central system ("dot") Hamiltonian i?systcm = 
J2a E'a\a){a\ expressed in terms of the system many-body 
eigenstates \a) and eigenenergies Ea, standard nonin- 
teracting free-electron Hamiltonian for the two mutu- 
ally biased (/x^ = eV/2, = —eV/2) leads i?ieads = 
'^k a=L R -^kacl^Cka, and the tunneling Hamiltonian 
Ht ^ T.ko.,ab [Taka)clja){b\+TUka)\b){a\cko]. 

Using the standard prescription^ we can write down 
the extended Liouville-von-Neumann equation for the 
generalized (non-Hermitian) density matrix p{t, x) — 

{xltXr) involving the counting fields Xa dual to the 
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number of charges passed through junction a (left or 
right) during time t. This equation reads (we set h = \ 
and e = 1 throughout the rest of the paper) 

ij^p{t, x) - H+{x)p{t, x) - pit, x)H-{x), (1) 

where the ^'dependent Hamiltonian is gener- 
ated by modifying its tunnel part as follows: 

H^ix) = ifsystem+i/lcads + Efea.a. K X)4a I «) (^1 + 

ri^(fc",x)|&)(a|cfca], withi3:T±(fca,x) ^ e^^^-^^T^^ika) 
and T±(/ca,x) = eT^x./27^,^(fca). 

For states of the whole system plus leads, we choose 
a basis of tensor products \ag) = \a) ® \g), with \g) 
a many particle state of the leads. We first eval- 
uate the time-evolution of system matrix elements 
Wab{t,x) = J2g{0'9\pit,x)\bg) and generalized system- 
lead coherences 4>ab{t,X\ka) = J2g{"'9\claPit,x)\bg) to 
which Wabit, x) directly couple. We then continue for the 

matrix elements of c^^cwa' p{t,x) and c\^c\,^,p{t,x)- 
At this level, we apply factorization and truncation 
conditions analogous to Ref. and, corresponding to the 
resonant tunneling approximation of Ref. 0, keep "non- 



diagonal matrix elements of the total density matrix up 
to the difference of one electron-hole pair in the leads" 
— i.e., for example, Y.g{o.g\c\^Ckacl,^,plt,x)\bg) ~ 
fkaJ2giO'9\cl,„,p{t,x)\bg) for ka ^ k'a', 

{ag\ckiaicl^a2^k3a3Pit,x)\bg) « 0, for kiai, ^202, and 
fcatta all different. We obtain a closed (though infinite) 
set of linear equations of motion for the reduced density 
matrix U'a6(i,x) and coherences 0a(,(i, X! ^Q^), which 
are more conveniently expressed in the Laplace picture 
{wab{z,x) = dte~^-*Wab{t,x) and analogously for 
(f>abiz, x; ka)). Assuming (l)abit = 0, x; ka) = and using 
the abbreviations = T^{ka,x), (t'ab = 4'abiz,x;ka), 
^ab = (l}li,iz,-x;ka), T'^k = T^f^{k'a',x) and anal- 
ogously for (j)'^^, (j)'^^, Tab, and T^^, our equations of 
motion read 

(iz ~ Ea + Eb)wab{z, x) - iwabit = 0, x) 

a' ,ka b' ,ka 

(2a) 



{iz - Ea+Eb+ Ek)(i>ab =Y^aa'fkaWa'b{z,x) -^'^b'bi^ ~ f ka)^ ab' {z , x) 



/ J aa 

a' ,k' a.' ^ka. 



La" [^a'-a'i^ ^ fk'a')4'a"b " a" fkaft^'ba"] + J2b' T f k' a' (p a' b' - T^'bi^ - fka)<P'b'a' 



iz — Ea' + Eb + Eka — Ek'a' 



y f+ - 

/ J a' a 

a' ,k' a'^ka. 

E Tb. 

h' ,k' a' ^ka 



Y.a" T'^,^„fk'a'<f>a"b -T^>a"fka4>a"b + J2b' T^' b'biX - fk' a')4>a'b' ~ T^^^i^ ~ fka)<i3'a' 



iz - Ea' +Eb + Eka + E, 



k' a' 



J2a' ^'a'a(l " fk' a')(l>a' b' " T+^' f kc4>'b' a' + Efc" T \, ^„ f k' a' (j) ab" " Ty,y{l - fka)(t)'b"a 



b' ,k' a' ^ka. 



iz — Ea + Eh' + Eka — Ek'a' 
Y.a' T''^^,fk'a'4'a'b'-T^^'fka(l>'a'b' +J2b" T' h"b'{^ - fk' a')4>ab" - Tf^'b,{^ - fka)'i>'ab" 



iz-Ea + Ey + Eka + Ek 



(2b) 



Eqs. pap and (j2b|) constitute the generalization of 
Eqs. (10) and (11) in Ref. i4j to case with the counting 
field; they are the main formal result of our paper and 
the starting point for the following studies. Analogously 
to the mean current,'^ charge conservation can be proven 
for all stationary cumulants using the method described 
m Ref. li!. Eq. ^ can be (at least formally) solved for 
4>ab in terms of Wab{z, x) and substituted into Eq. ((2a)) to 
give a closed non-Markovian generalized master equation 
for the reduced density matrix Wab{z,x) only. Its evolu- 
tion kernel VV(2, x) (s ee below) can then be used in the 
machinery of Refs.Tlsl andllGl to produce the current noise 
and higher-order cumulants of the full counting statistics. 



For ease of notation, we will only explicitly demonstrate 
this general procedure on the following example. 

Single-resonant-level model. — We consider an 
archetypical model of spin- less electrons and a single level 
forming the system TJsystom = ^^d|l)(l|. The only non- 
vanishing Tiia(ka) is rio(fca) = tka and consequently 



^ 10 



e=F^^°/%a, = e±*x./2t*^. Defining analo- 



10 — ^ '-ka 

gously to Ref. |4| (apart from the factor of 2tt) the quan- 
tity Ba{E) = 2TrJ2ktlaMz,X;ka)6{E ~ Eka) (with 



Ba{E)^B:{E) 



and introducing the conventional 



tunneling rates TaiE) = 27Tj2k\tka\'^^iE - Eka) and 
T{E) = TLiE) + Ti^iE), we get from Eqs. ^ and ^ 
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/dE — 

/fl E - 
— [Ba^iE) ~ B^E)] , 



E - Ed + iz - 



dE' T{E') 
E-E' + iz 



B^{E)^e-''^-/^fo.{E)V^{E) (^yjoo{z,x) +Y.e'''-'^' J 



dE' Ba,{E') 



2tt E' ~E-iz 

dE' Bo^'iE') 
'2^ E' - E-iz 



(3) 



We have verifiedi^ that, for x = 0; these equa- 
tions coincide with the RTA obtained by the real-time 
diagramatics/ 

From now on, we focus on a specific example of an 
infinite bias, i.e., Jl^E) = 1, fniE) = 0, with con- 
stant left tunnel rate TlIE) = Tl, and with the right 
tunnel rate having the Lorentzian energy dependence 

^r{e) 



2\n\' 



Sr 



This model has the ad- 



vantage of being analytically solvable for both the ex- 
act solution as well as the RTA equations, yet it ex- 
hibits nontrivial non-Markovian dynamics. This model 
is equivalen t^^i^^ to a noninteracting double dot with two 
electronic levels at Ed and Er , mutually coherently cou- 
pled by the transfer amplitude il, with the left level at 
Ed coupled by Fl to the left (filled) lead, and with the 
right level at Eji coupled to the right (empty) lead by an 
energy-independent rate 26 r. The double dot model can 
be described exactly by a Markovian generalized mas- 
ter equation, which can be then projected onto the left 
dot with the level Ed only resulting in an exact non- 
Markovian kernel VVoxact(2^, x)- Its explicit form is too 
lengthy and cumbersome to be presented here so we just 
use it for the exact evaluation of the noise and skewness 
shown in Fig. [TJ 

For the Lorentzian model, we can analytically solve the 
2vN equations ([3]) , which is in itself remarkable and does 
not work at finite bias. This solution is obtained from the 
ansatz BuiE) = hR{E)/{E - t - i^r), Bl{E) = h^iE), 
where bL^R{E) has no singularities in the upper complex- 
E half-plane. When this ansatz is inserted into the last of 
Eq. Q adapted to our example, one obtains, consistently 
with the assumed analytic structure in E, the following 
solution: 



Bl{E) = c-'^^/^Tl 
Br{E) = -Vr{E) 



E + iz- i:{E + iz) 



E- 



Y:{E + iz) 



(4) 



withe = ER-Ed and £(£') = -iV L/2+\n\'^ / {E-e + i5), 
the self-energy due to coupling to the leads corresponding 
to the sum of the tunneling rates. Using the relation 



I 

bR{e + i5R) = Yiu\E^e+isAE-{^ + i5R)]BR{E) stemming 
from the above ansatz for the determination of 6;j(e-|-i(5/j) 
from the second of Eq. Q and inserting the resultant 
form of Eq. ^ into the first two lines of Eq. ([3]), we 
finally arrive at the non-Markovian kernel in the 2vN 
approximation. 



W2vn(2,X) = 




- e^+l(z)['i(z) + \n\'^/(&n+z/2)\ 

2|np7(^) 

«"+7(z)l7(^)+z+|0| V(<5h+z/2)1 . 



(5) 

with 7(z) = 5r + Tl/2 + z. 

For zero counting fields, this kernel is identical to the 
exact one, W2-w^{z,Q) = VVoxact(2, 0), with clear conse- 
quence that the time-dependent occupations and station- 
ary current are also exact in accordance with previous 
findings.— 'ii'i^ The non-equilibrium rates are captured 
correctly including the effects of broadening of the level 
Ed due to the coherent coupling to the leads. However, 
already at first sight the counting field is accounted for 
in a rather primitive manner analogous to the lowest- 
order sequential tunneling model, i.e., the exact rates are 
just multiplied by the exponentials with counting fields 
(compare with Eq. (19) of Ref. [ish . For non-Markovian 
kernels, this prescription is potentially problematic and 
the exact kernel indeed contains the counting fields also 
in the denominators of the expressions for the rates ensur- 
ing consistency. Therefore, the 2vN kernel with counting 
fields is not exact and the noise and higher-order cumu- 
lants it yields are not correct. 

This is explicitly demonstrated in Fig. [U where we 
compare the exact solution for the noise and third cumu- 
lant with the 2vN/RTA and its Markovian limit obtained 
by using z in the kernel (O We also compare with 
the first-order von Neumann approach (IvN) equivalent 
to the standard lowest-order GME described in Sec. IIIA 
of Ref. [1^ and its Markovian limit. We fix the right tun- 
nel rate rfl(O) by adjusting while changing the param- 
eter 5r controlling non-Markovian behavior and observe 
the performance of various approximations for increasing 
degree of memory with decreasing 5r. The Markovian 
limit of the IvN is left constant by this procedure and 
serves as a reference, while all other solutions respond to 
the change of 5r. Obviously, none of the approximations 
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FIG. 1. (Color online) Second (top) and third (bottom) nor- 
malized cumulants of the current as functions of the width 
&B. of the Lorentzian right tunneling rate measuring the de- 
gree of non-Markovianity. Compared are the exact solution 
(solid black line) and four tested approximations: 2vN/RTA 
(long-dashed red line), its Markovian limit (short-dashed blue 
line), the first-order von Neumann/standard GME (long-dot- 
dashed green line), and its Markovian limit serving as the 
reference (short-dot-dashed magenta line). Parameters: 

e = rfl(o) = 4ri. 



reproduces the exact solution for strong-enough non- 
Markovian dynamics; all approximations perform quite 
badly for strong memory at (5^ ~ rfl(O) for both cumu- 
lants with extreme errors in the non-Markovian versions. 
These errors are typical for non-Markovian master equa- 
tions which have no Lindblad form and hence do not 
guarantee the conservation of probabilityi^. However, in 
the intermediate memory regime (S/j > rfl(O) the non- 
Markovian version of 2vN/RTA is clearly by far the best 
approach. 



Conclusion. We have extended the second-order von 
Neumann approach for GME kernels to the case with 
counting field which has enabled us to study its perfor- 
mance in the evaluation of the electronic current noise 
and higher cumulants. We have shown on an analytically 
solvable example of a single-level junction with struc- 
tured coupling to the leads that 2vN is not exact for the 
noise and higher-order cumulants and that it may fail 
significantly for highly non-Markovian dynamics. How- 
ever, for intermediate degrees of memory, 2vN seems to 
perform the best out of considered standard approxima- 
tions. Thus, when supplemented with efficient numerical 
implementation, it would be a method of choice for evalu- 
ation of noise and full counting statistics for systems with 
strong interplay between correlations and coherence. 
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